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ABSTRACT 

Context. In the past years, several observations of AGN and X-ray binaries have suggested the existence of a warm (T ~ 0.5 - 1 keV) 
and optically thick (tcor ~ 10 - 20) corona covering the inner parts of the accretion disk. These properties are directly derived from 
spectral fitting in UV to soft-X-rays using Comptonization models. However, whether such a medium can be both in radiative and 
hydrostatic equilibrium with an accretion disk is still uncertain. 

Aims. We investigate the properties of such warm, optically thick coronae and put constraints on their existence. 

Methods. We solve the radiative transfer equation for grey atmosphere analytically in a pure scattering medium, including local 
dissipation as an additional heating term in the warm corona. The temperature profile of the warm corona is calculated assuming it is 
cooled by Compton scattering, with the underlying dissipative disk providing photons to the corona. 

Results. Our analytic calculations show that a dissipative thick (Tcoi in the range 10-12) corona on the top of a standard accretion 
disk can reach temperatures of the order of 0.5-1 keV in its upper layers provided that the disk is passive. But, in absence of strong 
magnetic fields, the requirement of a Compton cooled corona in hydrostatic equilibrium in the vertical direction sets an upper limit on 
the Thomson optical depth Tcor 5 5. We show this value cannot be exceeded independently of the accretion disk parameters. However, 
magnetic pressure can extend this result to larger optical depths. Namely, a dissipative corona might have an optical depth up to ~ 20 
when the magnetic pressure is 100 times higher that the gas pressure. 

Conclusions. The observation of warm coronae with Thomson depth larger than =: 5 puts tights constraints on the physics of the 
accretion disk/corona systems and requires either strong magnetic fields or vertical outflows to stabilize the system. 
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1. Introduction 

Many successful models of the broad band spectra of accreting 
black holes contain a contribution from a moderately hot, opti¬ 
cally thick layer on the top of the relatively colder accretion disk. 
Such a layer is frequently postulated when the broad band spec¬ 
trum of an accreting system resembles a power law in the soft 
X-ray band, with the photon index above 2. 

The standard theory of the Shakura & Sunya^ l 1973 here¬ 
after SS disk) accretion disk model predicts that the spectrum 
of an accretion disk is well modelled as a multi-color black 
body, showing an exponential cut-off at high frequencies con¬ 
nected with the maximum of the temperature in the disk atmo¬ 
sphere. Some exceptional quasars/active galaxies and specific 
spectral states of X-ray binaries selected for the determination 
of the black hole spin show such a thermal cut-off ( [Steiner et al. 
201 Ij l, but in many cases the spectrum continues as a relatively 
steep/soft power law. 

The hard (2-100 keV) X-ray spectra of radio quiet AGN is 
generally characterised by a flat power law shape sometimes cut¬ 
off around 100 keV ( Jourdain et al.|1992[|Maisack et al.|1993[ 
Perola et al. [2002 1 Ballantyne 2014 Malizia et al. 2014 and 


references therein) and the presence of reflection components 
(iron line, reflection hump) is commonly observed ([Pounds et al.j 


19901. The soft (below 2 keV) band is generally characterised 


by an excess with respect to the extrapolation of this hard X-ray 
power law. This is the so-called soft X-ray excess. When fitted 
with a power law, it shows a steep (T > 2) spectral shape. The 
origin of this component is still unknown. It could be equally 
fitted by a blurred ionised reflection ( Crummy et al.[[2()06[ l, a 
blurred ionised absorption ( Gierlinski & Done|2004[ l or by ther¬ 
mal Comptonized emission in a moderately hot, optically thick 
layer possibly located on the top of the colder accretion disk 


([Walter & Fink|1993[[Magdziarz et al.|1998[ Done et al.|2012 

[Petrucci et al.|2()13[ l. 

The quasar composite spectra are well explained by such 
component, with the photon index T ~ 2.5 ([Laor et ah] 1997[ 


Elvis et al. 2012 1 . The whole class of Narrow Line Seyfert 1 
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galaxies has also similar soft X-ray slopes. Such specific spec¬ 
tral element is observed in galactic sources in the Very High 
State and Intermediate State (e.g. [Gierlinski & Done|20()3] l. The 
roughly power law shape of this component as well as the corre¬ 
lation observed between the UV and soft-X-ray bands suggests 
Comptonization as a mechanism responsible for this emission. 
Furthermore, observations show that the significant fraction of 
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the bolometric luminosity of the accreting system (~ 30 - 50%) 
is carried out by this emitting layer ( Vasudevan et al.|2014[ ). 

Usually, it is postulated that the hot medium responsible for 
this radiation forms a skin or a corona above the inner parts of 
the accretion disk. This skin or corona is optically thick in many 
models of specific objects, with Thomson optical depth of the 
order of 2 - 20. The observed slope of the soft X-ray spectrum 
does not determine the optical depth of the scattering medium 
since it depends on the Compton y parameter, i.e. a combination 
of the optical depth and the temperature. However, the fact that 
an unscattered disk component is not required in the fit or the 
direct detection of the turn-off imply rather low temperature and 
high optical depth. 

Such solutions are usually discussed in the context of spe- 
cihc observational data. White & Holt ( 1982) l required an opti¬ 
cally thick corona to explain the temporal and spectral properties 
of neutron star sources 4U 1822-37, 4U 2729H-47 and Cyg X-3 
(see also Bayless et al.|2010 1 , although recent papers (e.g. laria 


et al.||2013|l argue that the corona is very optically thin and the 


direct view to the neutron star in 4U 1822-37 is blocked by the 
outer rim of the accretion disk. 

Magdziarz et al. ( 1998|l postulated the presence of ; 


optically thick Comptonizing medium to ht the UV-Soft X-ray 
spectrum of NGC 5548. However, they did not expect it to co¬ 
exist with the cold disk as a vertical layer; instead, they sug¬ 
gested this medium as a radial transition region between the cold 
outer disk and a hot inner flow. Zhang et al. ( 2000| l modelled the 
spectrum of GRO J1655-40 with a warm layer at T=1.0 keV 
with the optical thickness t = 10 located above a cold accre¬ 
tion disk. Indeed, they have considered three vertical layers; a 
cold disk, a warm skin and a hot corona. [Janiuk et al.| ( |2001| l 
required the presence of a warm corona of optical depth equal 
to 12 in order to ht t he soft X-ray spec trum of a quasar/NLS 1 
object PG 121 1 h- 143. Zycki et al. (20011 found coronal temper¬ 
ature ~ 5 keV and optical depth ~ 3 from their hybrid models 
of soft states of X-ray binaries; GS 112468 and GS 2000-1-25, 
and higher values of optical depth, ~ 10, were implied for some 
of the Very High State data sets. Petrucci et al. ( 2013[ ) for the 
Seyfert 1 galaxy Mrk 509, found soft corona with temperature 
~ 0.5 keV, and large optical depth ~ 20. 

Monte Carlo models of hot coronae with large optical depth 
were calculated by Czerny et al. (20031 in order to explain the 
broad-band spectra of quasar composite spectra. Ton SI80, Mrk 
359 and PG 1211-1-143, i.e. high accretion rate AGN. [Kubota| 
& Done ( 2004| l obtained good spectral hts with the temperature 
of the order of 10 keV and optical depth ~ 2 in the two data 
sets for VHS of XTE J1 550-564. IJmerar] successfully 

htted the broad band spectra of 51 AGN with a model of ac¬ 
cretion disk thermal emission, a low temperature optically thick 
Comptonization and a hot optically thin corona. In their hts the 
electron temperature in the thick corona was in the range of 0.1 
- 2 keV, and the optical depth from 4 to 40 in different objects. 
The model of the optically thick, low temperature corona sur¬ 
rounding the cold disk was also successfully applied to model 
the spectrum of a ULX source IC 342 ( |Ebisawa et al.|2003 i al¬ 
though the model was considered rather unphysical by the au¬ 
thors. 

Postulating an optically thick hotter medium sandwiching 
the colder disk inside is in apparent conhict with the results for 
the radiative transfer in the diffusion approximation. A tempera¬ 
ture inversion is usually only obtained in the optically thin zone, 
whereas the temperature rises towards the interior of the celestial 
body in the optically thick zone. The best example is the solar 
corona. Thus the question arises of whether the cold disk would 
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Fig. 1. Sketch of the dissipative slab corona atop the non- 
dissipative atmosphere of the accretion disk. 


not heat up to reach the same temperature as the optically thick 
part of the corona? 

In this paper, we address this question using a very simple 
analytical model. We consider the vertical structure of the accre¬ 
tion disk/corona that is sketched in Pig.[T] and investigate the ra¬ 
diative and pressure equilibrium of the upper layers of the accre¬ 
tion flow as a function of the fraction of the accretion power that 
is dissipated in the corona. We show, that the disk embedded in 
the hotter optically thick medium indeed does not heat up much 
more than in the case of an optically thin surrounding corona 
discussed in a basic paper of Haardt & Maraschi] ( |1993 i. We 
check conditions for which an optically thick, Compton cooled 
zone can exist in hydrostatic equilibrium with a specihed under¬ 
lying colder disk. The paper is organized as follows; in Sec. 
we describe the solution of radiative transfer equation for the 
grey atmosphere with additional heating, in Sec. |^the tempera¬ 
ture structure is presented, in Sec.j^the hydrostatic equilibrium 
is solved, and in Sec. we show when corona is dominated by 
Compton cooling. All results are discussed in Sec. 


2. Grey optically thick scattering medium with 
dissipation 

We are interested in the case of a grey atmosphere with pure 
scattering (i.e. we neglect emission and absorption). We assume 
that the atmosphere is optically thick, and therefore adopt the 
Eddington approximation in the whole medium. We allow for an 
additional dissipation which heats the corona. Eor simplicity, we 
assume that the input energy rate per unit optical depth (and solid 
angle), Q is uniform from the surface (t = 0) to the base of the 
corona (r = Tcor)- This corona is located above a cold accretion 
disk which itself dissipates energy due to the accretion flow. The 
assumption that all dissipated energy is emitted in the form of 
radiation allows us to derive a simple fully analytical solution 
for the local vertical temperature structure. At the surface of the 
corona the escaping radiation flux E™ is the sum of the flux 
generated through internal dissipation inside the accretion disk, 
7^di.sk5 and that produced through dissipation in the corona, 

= ^disk + . (1) 

The assumption that the heating of the corona is uniform gives; 

Ecor = 47rgrcor ■ ( 2 ) 


Additionally, we define a convenient parameter; 


A ^ 


(3) 
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which is the fraction of the total accretion power that is dissi¬ 
pated in the corona. 

The frequency-integrated radiation transfer equation with an 
additional energy input Q can be written as; 

dl 

(4) 

dT 

where ji is the cosine of azimuthal angle. The optical depth, t 
is measured downward, from the top of the corona toward the 
disk, /(/r, t) is the radiation intensity, and J{t) is the mean inten¬ 
sity. The third term on the right hand side modihes the source 
function, S , and describes the increase in the photon energy due 
to the dissipation within the corona (S - J + Q). Following the 
standard Eddington approach, we can derive the solution of the 
radiative transfer equation by calculating its hrst moments, i.e. 
integrating over solid angles. The zeroth moment gives: 

H{r) = -tQ + Cl, (5) 



The integration constant Ci represents the Eddington flux at the 
top of the corona (r = 0 ), which is, by definition, Ci = E‘°‘j,/ 47 r. 
Using Eqs.[^and[^ the Eddington flux profile can then be written 
as; 


rtot 

H(t) = 

An 



(6) 


We note that at t = Tcor where the corona touches the cold 
disk, the downward flux corresponding to the illumination of 
the disk by the corona cancels out the upward flux of repro¬ 
cessed/reflected radiation from the disk. Therefore the net ra¬ 
diation flux at Tcor is only that caused by internal dissipation in 
the cold disk Edisk- 

The first moment of the radiation transfer is; 


K{t)^ 




+ C 2 . 


(7) 


Fig. 2. Typical radiation structure of the accretion disk corona 
and upper disk atmosphere. The upper panel shows the assumed 
dissipation profile, the middle and lower panel show the resulting 
radiation flux and pressure profiles, respectively. 


The solution of Eq. with Q - 0, implies a constant 
Eddington flux below the corona. Then the value of the flux at 
Tcor sets the value of H everywhere in the disk atmosphere; 

H(t > Tcor) = ^ = ^(1 -X) . (9) 

4;7r 47r 

The first moment of the radiation transfer then gives; 

plot 

^:(T>Teor)= ^(1 -T)t+C3. (10) 


In the Eddington approximation, we accept K - 7/3 at every 
optical depth across the medium. In addition, at the corona sur¬ 
face we have only outgoing radiation flux, as in standard stellar 
atmosphere, so we have the condition 7(0) = 2H{Q) which al¬ 
lows to determine the constant C 2 . Setting t = 0 in Eq.j^an d 0 
we get C 2 = E'cc/hTT, so that the mean intensity as a function 
of the optical depth in the optically thick corona is given by the 
expression: 


7(t) = 


3F'. 


An 


2 

n + T- 


XT 

O'T 

^ ‘ rn 


( 8 ) 


Equations and 1 ^ are valid only in the warm corona (i.e. 
for T < Tcor) and they are largely independent of the underlying 
accretion disk structure. However the same formalism can be 
used to extend these solutions deeper in the disk atmosphere in 
order to investigate the effects of the presence of the corona on 
the upper layers of the accretion disk. 

As in the standard SS 1-D disk, all flux is generated close 


to the equatorial plane (Rozanska 


19991, we can neglect dissi¬ 


pation in the disk atmosphere, and assume that all the disk flux 
is generated below this layer at deeper optical depth (t > 10^), 
see Eig. [T] Here we consider only the properties of the upper 
non-dissipative atmosphere. Eor simplicity, we will not solve the 
flux dissipation equation down to the midplane as this would not 
change our final result of temperature profile on the border of the 
warm corona and the disk. 


The constant C 3 is determined from the condition of continuity 
of K (and 7) at r^or: 


C 3 


rtot 

acc 


47r 



( 11 ) 


And Anally the mean intensity in the disk atmosphere is; 


T rtot 

77 ^ \ acc 

> Tcor) = —;- 


X 2 

(1 -T)t+ + 3 


( 12 ) 


We note that in absence of corona (r = 0) this equation reduces 
to the standard mean intensity profile of grey atmospheres: 


, _ 3Edisk ( 

■^d.sk - ^ ^ 


2 

T -H - 

3 


(13) 


In the presence of a warm corona, there is an additional com¬ 
ponent 7cor to the mean intensity of the disk that is due to the 
illumination of the atmosphere by the corona, 7 = 7disk + Jcor, 
which is: 



(14) 


The typical vertical structure of the radiation properties of the 
corona/disk atmosphere system is sketched in Eig. 
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Fig. 3. Temperature profiles in the disk/corona system for the 
various values of the parameters x, fcor of the six fiducial mod¬ 
els reported in Table Each curve is labelled with its model 
reference number, corresponding to that given in Table [T] 


3. The temperature profile 

The next step is to determine the temperature prohle across the 
disk and corona using mean intensity field. The temperature pro¬ 
hle is derived assuming that matter is in equilibrium with the ra¬ 
diation held. The resulting temperature depends on the radiative 
cooling mechanism operating in the gas. 


3.1. Temperature profile in the warm corona 

The radiative transfer solution for grey atmosphere does not 
specify the temperature prohle in a purely scattering medium 
which is the case of the soft corona, i.e. for t < Tcor- However, 
we can obtain the temperature prohle taking into account that the 
scattered photons must cool the corona through Comptonization, 
in order to have the thermal heating/cooling balance. 


Table 1. Parameters of the hducial solutions illustrated in Figs.j^ 
and|5 For all these models, the total accretion hux is set to = 
3.3x10*^ erg s“' cm“^. The values of Tcor were choosen to be the 


maximum possible Thomson depth of a Compton cooled corona 
in hydrostatic equilibrium for the corresponding JS^, Q and;^f (see 
Sect. 1^. Tav is the resulting average temperature of the corona 


estimated using Eq. 17 


Model number 

Pm 

9 

A 

^cor 

kT,, (keV) 

1 

0 

0 

0.98 

5.21 

3.91 

2 

50 

0 

0.98 

19.9 

0.42 

3 

50 

5 

0.98 

8.76 

1.68 

4 

50 

0 

0.4 

15.9 

0.23 

5 

50 

2 

0.4 

7.08 

0.88 

6 

50 

0 

0.02 

9.28 

2.68 xlO-2 


When cooling is dominated by inverse Compton scattering in 
the Thomson regime by a thermal population of sub-relativistic 
temperature, the balance between cooling and heating reads; 


J(t) 


4^7’cor(T) 


ms-c^ 


Q 


(15) 


where Tear is the corona temperature, c - the velocity of light, k 
- the Boltzmann constant, and - the electron rest mass. 

Since J is an increasing function of the optical depth (see 
Eq. and Q is assumed to be a constant, the electron tempera¬ 
ture decreases with t, showing a temperature inversion; 


^Tcoi(t) = 


Xm^ 

l^Tcor 



(16) 


The temperature remains hot from the coronal surface down 
to T = 2/3, then it slowly decreases. Depending on the amount 
of dissipated energy, the soft corona can be quite hot down to 
moderate optical depth. The vertical temperature profile of the 
soft corona above the cold disk is determined by only two pa¬ 
rameters; the total optical depth Tcor of the skin/corona, and the 
fraction x of flux dissipated inside the corona to the total flux 
dissipated inside the disk and corona. In Fig. we show exam¬ 
ples of full temperature profiles for different values of these pa¬ 
rameters. Note, that this profile strictly comes from the radiation 
properties and does not imply that the hydrostatic equilibrium is 
satisfied in this multi-zone structure. We address this issue in the 
section below. 

One-zone Comptonization models such as those used to de¬ 
rive the observational properties on the warm corona only con¬ 
strain the average temperature of the skin, without any constrains 
on the temperature stratification. Since, in our models tempera¬ 
ture changes with optical depth, we should compare the value of 
the average temperature of the warm skin weighted by its optical 
depth with that determined from observations; 



Tcoi(T)dT = 


XmeC^ 

12^MTcor 


1 + 

u~l I 

1 - 

^ M+l ' 


where 


(17) 


u - 



(18) 


In the limit of large Tcor, the average temperature given by Eq.fTT] 
can be approximated within 10 percent (for Tcor > 3) as; 


- 


12 tL 


In 


3 tc, 


-T 


(19) 


The dependence of the average coronal temperature on Tcor is 
plotted in Fig. for several values of the x parameter. We can 
clearly see, that a warm skin, cooled by Comptonization, can 
be produced even for large values of the coronal optical depth 
provided that most of the accretion power is dissipated in the 
warm corona. The coronal temperature however decreases with 
decreasingand increasing Tcor- 

As a consequence, at small i.e. strong disk dissipation, a 
hot corona with pure scattering is only consistent with the most 
moderate observations (i.e. Tcor of the order of a few). From 
the energy equilibrium requirement we can produce a layer of 
Tcor ~ 0.5 - 1 keV and large Thomson depth Tcor <15 only if 
most of the accretion power is dissipated in the layer. This ap¬ 
pears to be consistent with the observational results of |PetruccI 
|et al.|(|2013|l in the case of Mrk 509. These authors infer such 
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parameters for the soft corona and also argue that the observed 
relative luminosity of the disk and soft corona implies that the 
disk is passive (i.e. x - The case of a passive disk provides 
the maximum achievable temperature for a given coronal depth. 
We note that some numerical simulations of accretion disks also 


show a stronger dissipation in the outer layers of the disk (Hirose 
|& Turner|2011| l. 

We also note that the method presented here is reason¬ 
ably accurate. Our results can be compared for example to the 
Monte-Carlo simulations presented in Malzac, Beloborodov & 
Poutanen (2001), in the case of a slab corona with Thomson 
depth of Tcor = 3 above a passive disk. They obtained an av¬ 
erage temperature ^ 9 keV (see their hgure 2) which is in ex¬ 
cellent agreement with the present results. On the other hand 
for optically thin, coronae, the temperature becomes mildly rel¬ 
ativistic, our approximations break down and the simple analyt¬ 
ical model underestimates radiation cooling. For Tcor = 0.5 and 
X - Eq. 17 gives kT^y ^103 keV while Malzac, Beloborodov 
& Poutanen (2001) obtain kT^y ^ 70 keV with their detailed 
calculation. 


3.2. Temperature profile in the disk atmosphere 

Deep in the optically thick atmosphere of the disk (t > Tcor) 
we can assume that radiation is fully thermalised {J - 
where cr is the Stefan constant). Using Eq. [T^ we obtain the 
following temperature structure: 


T rtot 
y4 _ acc 


n \ ^ X'^cot 


This expression can be rewritten as: 


^ atm “ ^ ^ disk 


2 

T + - 

3 


71J c( 


( 20 ) 


( 21 ) 


The hrst term on the right hand side of Eq. [^corresponds to the 
standard temperature structure for fully thermalised grey atmo¬ 
sphere in the Eddington approximation. Tdisk is the disk effective 
temperature in absence of warm corona, calculated from the in¬ 
trinsic disk flux Fdisk. The constant second term represents the 
increase in disk temperature due to the coronal illumination. 

Unlike Tcor, the temperature prohle of the disk atmosphere, 
Tatm, depends on the accretion flux F'cc. Using standard ac¬ 
cretion disk theory, Face estimated as a function of the 

mass of the black hole, Mbh ,and Eddington luminosity fraction, 
m - L/Ledd. a given radius R - rGM^nlc^' 


F“c-8x102'> -4 

tvj yJ 


erg s ’ cm ^ , 


( 22 ) 


where m - Mbh/Ms, / = 2 r;(l- yrifr) and r, is the inner radius 
of the disk expressed in gravitational radii. 

In Eig. [^ the prohle in the disk atmospher e was calc ulated 
for a black hole of mass Mbh = 1-4 x IO^Mq ( Liu|l983 i, at 10 
gravitational radii from the black hole, and at an accretion rate 
equal to 2% of the Eddington accretion rate, and we set r, = 5. 
Eor these parameters the accretion hux is F'ec = 3.3 x 10*^ erg 
s ' cm“^. As can be seen on Eig.[^ the temperature prohle in the 
upper layers of the disk battens and departs from the standard 
grey atmosphere temperature prohle due to the strong coronal 
illumination only when most of the power is dissipated in the 
warm corona. 

The absorption of coronal photons by the cold disk only 
slightly increases disk effective temperature which remains sig- 
nihcantly lower than that of the corona. Due to additional dissi¬ 
pation in the soft corona, all temperature profiles show a strong 



Fig. 4. Average coronal temperature vs. optical depth. The full 
curves show the dependence of the average coronal tempera¬ 
ture on the total optical depth of the corona for various values 
of the fraction of total power dissipated in the corona, x^ la¬ 
belled. The dashed curves represent the hydrostatic equilibrium 
solutions providing the largest possible Thomson depth of the 
Compton cooled corona for a given x and magnetic to gas pres¬ 
sure ratio jSni (see Seetj^. Each of the dashed curve shows the 
track of the solutions in the kTav-Tcor plane for a hxed value of 
jStn (as labelled) when x is varying. 


temperature inversion in the disk/corona system. At the transi¬ 
tion between the disk and corona there is a discontinuity due to 
the change in cooling mechanism. The amplitude of the temper¬ 
ature jump can be estimated as: 


T Sim 


2 / 
"H ^cor 1 


Tcor 

T=r„, \mGl X 

3 \ 

2 / 


This ratio remains lower than unity over a very broad range of 
black hole masses, mass accretion rates and disk radii. In prac¬ 
tice, only when x vanishes, the temperature of the disk at Tcor 
can become comparable or even hotter than that of the corona. 

We note that in reality the change in cooling mechanism 
might not be as brutal as we have assumed here and the temper¬ 
ature discontinuity could be smoothed. The computations of the 
disk/corona transitions are difficult but it appears however that 
the temperature drop is always very sharp even if all the radia¬ 
tion processes are fully taken into account (jMadej & RozahsE 
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2000bt Ballantyne et al.|[2001t [Nayakshin & Kallman|pOOT| 
Rozahs ~a et al.|2002|l. 


4. Hydrostatic equilibrium for the corona/disk 
system. 

The question arises if such a two-zone system can he in hydro¬ 
static equilihrium. Here, we derive analytical formulae for the 
pressure profile in the disk-corona system. We use the standard 
equation of vertical hydrostatic equilihrium in geometrically thin 
disk ( |Madej & R6zanska||2000a I. The total pressure, P, is the 
sum of the gas pressure Pgas, the radiation pressure Prad, and the 
magnetic pressure Pmag- Locally, the total pressure has to bal¬ 
ance the gravitational force: 


dP„ 


dP„ 


dr 


dr 


1 GMbh flfPrad 

dT ’ 




(24) 


The first term in the right hand side of Eq. is the ratio of 
the gravitational to radiation pressure force at the surface of the 
corona. For a corona in hydrostatic equilibrium this ratio is nec¬ 
essarily larger than unity and consequently ^ > 0. The half of 
the disk thickness Zdisk is a crucial parameter which controls the 
hydrostatic equilibrium in the vertical direction ( [Rozanska et al. 
\999\ . The case Q - 0 gives the minimum disk thickness, for 
which the total disk pressure can be balanced by the gravitational 
force: 


y^min _ 
■^disk “ 


K F™ 
^es ^ acc^'^ 

GMbh c 


3 GMbh 
2 


fm. 


(31) 


For a thinner disk, the matter will be outflowing from the system 
( jWitt et al.|19^ . 

The density profile clearly depends on the assumed disk geo¬ 
metrical thickness Zdisk. and we can derive this density from the 
equation of state as: 


where G is the gravitational constant, and ACes - the Thomson 
scattering cross section. In a general approach the total opac¬ 
ity should be taken into account, but for simplicity, we take only 
Thomson scattering into account and set /Cgs = 0.34 cm^ g“'. 
In order to obtain analytical solutions, we will assume that the 
warm corona is geometrically thin compared to the scale-height 
of the disk Zdisk, so that the vertical distance to the equatorial 
plane, z , can be considered a constant z - Zdisk- This implies 
that the gravitational force is constant along the vertical direc¬ 
tion inside the corona and the upper layers of the disk. 

For the grey atmosphere, the radiation pressure gradient de¬ 
pends on the flux expressed in Eq.j^for the corona and Eq.j^for 
the disk atmosphere: 


dPxaA 

dT 



(25) 


P = 


F'«h 


kFcor (1 +y6m)c 




XT 

'2.Trn 


(32) 


where fiis a mean molecular weight assumed to be 0.5, and mu 
is the mass of hydrogen atom. The density increases with t and 
from Eq. 32 we see that setting the disk parameter ^ = 0 (or 
equivalently Zdisk = -^disk^ minimises the density in the corona. 
This is illustrated in Fig. which displays examples of density 
profiles obtained for the parameters listed in Table [T] 


4.2. Pressure and density profiles in the disk atmosphere 

The pressure and density profiles in the disk atmosphere below 
the corona (t > Xcor) can be estimated in a similar way. The 
radiation pressure profile is given by Eq. 12 


4.1. Pressure and density profile in the warm corona 


ptot 

Pr.d = — 

C 


(1 -A)T + ^Tcor 


2 

3 


(33) 


In the corona, the radiation pressure profile is obtained directly 
from Eq. 


47tJ 

3c 


plot 

acc 

C 




(26) 


We assume a uniform magnetic to gas pressure ratio /3ni' 


P = 

-* mag 


Bi 


Stt 




gas 


(27) 


Solving the equation of hydrostatic equilibrium ( |24| l, we find an 
expression for the gas pressure structure assuming, as a bound¬ 
ary condition, that Pgas(T = 0) = 0: 


P 


gas — 


ptot 

acc 


(1 +^m)c 




(28) 


where the constant ^ represents the ratio of the pressure forces 
of the gas and magnetic field to that of radiation at the surface of 
the corona: 


& - 


BP gas ^ dPjjxag \ ^ dP rad 
dr dT I dT 


(29) 


The radiation pressure force dominates the support of the corona 
at all depths for Q < \ — 2x. 0 can also be expressed as: 


GMbh cZdisk 


(30) 


The pressure equilibrium equation (24 1 is solved using Eq.j^and 
assuming continuity of pressure at the disk/corona transition: 


p _ acc 

“ (1 +ySJc 




X 


(34) 


The density profile follows: 


P = 


pm-a 


kTaXm (1 +y6m)c 


+X)t-^t a, 


(35) 


Fig. 1^ shows some examples of density profiles around the 
disk transitions for our fiducial value of the total accretion flux. 
In all cases, we observe a discontinuity in density at Tgor which 
has the same amplitude as the temperature jump discussed in 
the Sect. |3.2| Due to the temperature jump and the condition of 


pressure equilibrium at the disk/corona transition, the disk atmo¬ 
sphere tends to be much denser than the corona. 


4.3. Limitations 

In Sect.j^we have estimated the properties of a warm corona and 
disk atmosphere in pressure equilibrium. Our proposed treat¬ 
ment presents the advantage of being very simple. The draw¬ 
backs of this simplicity are some limitations that we now briefly 
discuss. 

First, althought the effects of the disk on the corona and disk 
atmosphere is taken into account via the 0 parameter, the present 
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Fig. 5. Density profile around the disk corona transition for the 
6 fiducial models detailed in Table [T] The curves are labelled by 
their model number. 


approach does not allow us to guarantee that there is indeed a 
disk solution below the atmosphere that both has the required 
Q parameter and connects smoothly to the atmosphere. A full 
calculation of the vertical stratihcation of the disk down to the 
mid-plane, including also dissipation, would be required in order 
to obtain such self-consistent solutions and calculate Zdisk from 
first principles. 

Also we note that for the density profiles with the lowest 
densities in the corona, our assumption of constant gravity may 
be inaccurate. Indeed for these prohles the scale height of the 
corona: 

^Tcor 

I {K^spT^dT, (36) 

^Tcor/2 

can be be comparable, or even larger than Zdisk- In this case the 
gravity at the surface of the corona is significantly larger than 
at the bottom. Taking into account the height dependent gravity 
would require a numerical resolution of the equilibrium which 
is out of the scope of this paper, but we can anticipate its effects. 
Indeed, if we assume that the disk/corona transition is located at 
a height Zdisk, the increased gravity force in the upper corona will 
necessitate a larger pressure in order to sustain the equilibrium, 
and as a consequence the coronal density will also be increased 
compared to our current estimates. In particular this effect may 
affect our results for small Zdisk (or small Q), which may under¬ 
estimate the pressure and and density in the corona by a factor 
of up to a few. 


Finally our calculation of the pressure equilibrium assumes 
that the opacity is dominated by electron scattering. If absorption 
becomes important, both pressure and density will be reduced 
compared to our simple estimates. We have checked a posteri¬ 
ori that for the fiducial models presented in Figs. and 1^ the 
Kramer free-free absorption opacity A-ff ^ 6 x 
g ' is negligible compared to /Ces both in the the corona an the 
atmosphere of the disk. For these models our estimated pressure 
profiles are not affected by the approximation of a pure scatter¬ 
ing medium. We stress however that since even in the disk atmo¬ 
sphere the medium is only weakly absorbing, the assumption of 
fully thermalised radiation used to infer the temperature profile 
of the atmosphere may break down close to the disk/corona tran¬ 
sition, making the transition much more gradual than our simpli¬ 
fied calculation suggest. 

A detailed investigation of all these issues is deferred to fu¬ 
ture works. 


5. Constraints from the requirement of a Compton 
cooied corona 


In the previous sections we have determined the temperature, 
pressure and density profile of the corona under the assumption 
that the dominant cooling mechanism is Compton scattering. 
This assumption was motivated by observational results based 
on the modelling the coronal emission with Comptonization 
models. We now have to determine the parameter regimes 
for which this assumption remains valid. Besides Compton 
cooling, the most efficient cooling mechanism is expected to 
be bremsstrahlung which must remain negligible compared to 
Compton cooling. Here we estimate the ratio of Compton cool¬ 
ing rate Ac = IbnkT j{m^c^)p {t) (in erg/s/cm^) to the 
bremsstrahlung cooling rate Ab = where B - 6.6 x 

10^° CGS units. This ratio must remain larger than unity across 
the soft corona. Using Eqs.[^and[T^we get the condition: 


Ac 

Ab 






2 \~^ 
x-^ \ 

2Tcqi- j 


> 1, (37) 


with the constant A — ^Ikm^ VT2Z ?p mn) - 57. The 

above ratio is a decreasing function of optical depth, there¬ 
fore the condition is verified in the whole corona if it veri¬ 
fied at T = Tcor- The maximum possible Thomson depth for a 
Compton cooled corona is obtained by setting the condition that 
the depth Tcor corresponds to the transition between Compton 
dominated to bremsstrahlung dominated regions, i.e. by solving 
Ac/AbIt=t„, = 1 - Ths coronal optical depth of the fiducial mod¬ 
els of Table leading to the profiles presented in Fig. and 
were determined in this way and correspond to the deepest pos¬ 
sible Compton corona for the given set of^, and Q. 

From Eq. we also see that the case of ^ = 0, or equiv¬ 
alently Zdisk = -^^sk’ gives the most favorable condition to have 
a purely Compton cooled corona because it corresponds to the 
minimum of the gas pressure and density. Therefore, in order to 
find regimes where our assumptions fail, it is enough to calculate 
this ratio at the base of corona for 0-0: 


^ ^ V8A(1 +p^) 

Ab T=r,„ rlor V2/A - 1 


(38) 


We stress that this equation represents the necessary condition 
for the Compton cooling dominance and does not depend on any 
disk parameters. 
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Fig. 6. The importance of Compton scattering over the 
bremsstrahlung at the base of corona as a function of Tcor 
(Eq. [3^ . Solid lines in each package are computed for various 
values of;^f = 0.18,0.26,0.40,0.57,0.67 and 0.86. Each pack¬ 
age is calculated for three different values of magnetic pressure: 
ySm = 0,10 and 100. The horizontal dashed line represents the 
case where Ac/Ab = 1 , while vertical dotted lines mark values 
of coronal optical depth for which it happens. 


Fig. 0 shows the Compton to bremsstrahlung cooling ratios 
as a function of coronal optical depth for three different values 
of magnetic pressure in the case 0-0. Eor each value of 
we consider various values of x to show the marginal influence 
of this parameter. We see that the maximum possible Thomson 
depth of the corona in the case 0 = 0 provides an absolute upper 
limit on the depth of the corona. This upper limit can be esti¬ 
mated as: 

- 5.4 (1 . (39) 


The corresponding average temperature of this ‘deepest’ possi¬ 


ble corona can then be estimated using Eq. 19 


kT,. 4 


xH2-xy^ 

(1 +y6m)5 


1 - 1.2 X 10^'* In 


x-H2-xy 


(1 +jem)2 


keV. (40) 


Eq. shows that when magnetic pressure is zero, 
bremsstrahlung becomes the dominant emission process as soon 
as Tcor ^ 5. Such solutions are therefore inconsistent with an un¬ 
magnetised, static, Compton dominated, hot, and optically thick 
corona. Relaxing one of these constraints might produce con¬ 
sistent solution. Eor instance, the corona might not be in static 
equilibrium. The case of outflowing coronae ( |Witt et al.|1997| l in 
the frame of our model will be investigated in a future work. 

Alternatively, additional magnetic pressure helps the gas 
pressure to balance gravity. Hence it produces solutions with 
lower density, i.e. with lower bremsstrahlung cooling rate. Eig. ^ 
maps the iso-contours of x and /?„, for the deepest possible 
corona, in the Tcor vs kTav plane. This diagram allows one to 
estimate the minimum value of the magnetic field pressure ra¬ 
tio required in order to produce a corona with a given optical 
depth and temperature. Eor instance, the warm corona recently 
observed in Mrk 509 by Petrucci et al. (2013), with optical thick¬ 
ness Tcor ~ 15 and temperature - 0.5 keV, implies magnetic 
to gas pressure ratio > 30. Eor Mbh = 1.4 x IO^Mq, a dis¬ 
tance of 10 gravitational radii, and an accretion rate equals to 


2% of the Eddington accretion rate, this implies that the average 
magnetic field in the corona must be larger than 10^ G. 

Since the calculations are based on the assumption of ^ = 0, 
they provide the necessary conditions for the Compton cool¬ 
ing condition to be satisfied. If Compton cooling dominance 
is not obtained, no change of Zdisk will save the situation, and 
if Compton cooling dominance is obtained, there is always a 
range of Zdisk > -^disk’ which the hydrostatic equilibrium is 
sustained, and density is low enough for Compton cooling to 
dominate over bremsstrahlung cooling. We note however that in 
the case of minimal density in the corona (0 = 0) our approx¬ 
imation of constant gravity force in the vertical direction can 
break down (see Sec. |4.1| l. As a consequence our simple calcula¬ 
tions may underestimate the density and bremsstrahlung cooling 
rate. Relaxing the assumption of constant gravity may reduce 
the maximum possible Thomson depth of the corona and/or re¬ 
quire even stronger magnetic fields to maintain the dominance 
of Compton cooling. 


6. Discussion and Conciusion 


In this paper, we put constraints on the existence of a warm, 
dissipating, optically thick, and Compton cooled corona in hy¬ 
drostatic equilibrium with a cold accretion disk. We neglect syn¬ 
chrotron photons, thermal conduction and ionisation in the warm 
skin, but we have checked that those processes are not important 
comparing to Comptonization. In our computations of the hy¬ 
drostatic equilibrium in the vertical direction, the radiation pres¬ 
sure component is fully taken into account, contrary to numerical 
simulations by |Schnittman et aTT (20131; Uzdensky ( 2013| l. 

Our simple analytical solution for the warm and dissipative 
corona above the cold disk shows that a stable temperature inver¬ 
sion is possible in the optically thick case, contrary to the intu¬ 
itive expectations based on the diffusion approximation. This is 
due to the fact that the solution of purely scattering atmosphere 
does not specify the temperature. It only gives the radiation den¬ 
sity, which rises with the optical depth, as expected. The tem¬ 
perature in the warm corona is determined a posteriori, from 
Compton cooling balance, and the temperature increases toward 
the warm skin surface. We have shown that such corona can 
reach temperatures of 0.5-1 keV for assumed values of constant 
dissipation in the skin of moderate optical depth (Tcor < 10). The 
most extreme parameters e.g. a coronal temperature 0.5 keV at 
optical depth 15-20, which is observed in the case of Mrk 509 
( Petrucci et al.|20l3] l, can also reproduced provided that the disk 
is passive (i.e. almost all of the accretion power is dissipated in 
the corona). 

Nevertheless, if this zone is in hydrostatic equilibrium with 
the cold accretion disk, the maximum optical depth of the 
corona cannot exceed ~ 5 without any additional magnetic 
pressure.This upper limit is independent of the disk parameter. 
Higher optical depth of the warm skin is possible if the gas pres¬ 
sure is lowered, by magnetic pressure, or possibly by mass out¬ 
flow. In this paper, we illustrate only the first case i.e. non zero 
value of magnetic field strength. When, the ratio of magnetic 
pressure to the gas pressure is 100, the maximal optical depth 
of the warm corona is around 20, which is consistent with some 
observations. 

We conclude, that in the absence of magnetic pressure, addi¬ 
tional dissipation in the outer layer is able to heat up the corona, 
but the requirement of hydrostatic balance with the disk, puts 
strong limit on the coronal optical thickness. This limit is inde¬ 
pendent on the accretion disk parameters, i.e. its accretion rate 
and the mass of the central black hole. 
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In this context, the X-ray observations of optically thick 
('rcoi> 5), warm coronae have a strong implication for the 
disk/corona system: either strong magnetic helds or vertical out¬ 
flows are required to stabilise the system. What is more, the 
simple conditions discussed in this paper are the minimum re¬ 
quirements for the existence of the thick corona, and further 
modelling of the disk/corona interaction may likely impose even 
more stringent constraints for the existence of such medium. 
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